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A systematic analysis of the moments of the fragment size distribution has been carried 

out for the multifragmentation (MF) of 1A GeV Au, La, and Kr on carbon. The breakup of 

Au and La is consistent with a continuous thermal phase transition. The data indicate that 

the excitation energy per nucleon and isotopic temperature at the critical point decrease 

with increasing system size. This trend is attributed primarily to the increasing Coulomb 

energy with finite size effects playing a smaller role. 



PACS number: 25.70Pq; 05.70. Jk 

The EOS collaboration has recently studied the multifragmentation (MF) of 1A GeV 
Au on carbon [|IHT0H- One of the important results was the possible observation of critical 



behavior and the extraction of associated critical exponents ]T],[3],|5|] . The values of these 
exponents were very close to those of ordinary fluids indicating that MF may arise from a 
continuous phase transition and may belong to the same universality class as ordinary fluids. 
Another important result was the successful description of the EOS MF data by statistical 
thermodynamical models, which describe quantum mechanically the MF of a charged nucleus 
|T0|-p!5[. In this paper we analyze the recent results for MF of lAGeV La and Kr on C M 
along with those previously reported for Au |],|3]|| in the manner proposed by Campi |16HT8" 



Our analysis provides the first experimental evidence for the evolution of the MF mechanism 
with increasing projectile size and for the effects of Coulomb energy and finite size. 

The reverse kinematics experiments and the analysis by which the equilibrated remnant, 
which undergoes MF, was separated from promptly emitted particles and the details of the 
determination of the remnant mass and excitation energy are given in our earlier publications 

ill- 

Campi [16-18] and Bauer 1 1^,20] first suggested that the methods used in percolation 



studies may be applied to MF data. In percolation theory the moments of the cluster 



distribution contain the signature of critical behavior [21]. The method of moments analysis 



was used by several groups [22-25] to search for evidence of the liquid-gas phase transition 
in MF. Thus for each event, we determine the total multiplicity of charged fragments, m, 
and the number of charged fragments n^, of nuclear charge Z and mass A 0. We calculate 
the k moments of the cluster size distribution given by 

M k (m) = ^ASH (1) 

where the sum runs over all masses in the event including neutrons except for the heaviest 
fragment. This quantity was instrumental in extracting the critical exponents in Au+C data 



|Ij]. It has been argued that there should be an enhancement in the critical region of the 
moments, M&, for k > r — 1, with critical exponent r > 2 [T^|T^|. For example, the reduced 
variance 72, i.e. the combination of moments given by 

72 = M 2 M /Ml (2) 

has a peak value of 2 for a pure exponential distribution, ua ~ e~ aA , regardless of the value 
of a, but 72 > 2 for a power law distribution, ua ~ ^4~ T , provided the system is large enough. 
Here Mi and Mi are the first and second moments of the mass distribution in an event and 
Mq is the total multiplicity including neutrons. 
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FIG. 1. 72 as a function of multiplicity from Au, La, and Kr systems 



We have calculated 72 event by event as a function of total charged particle multiplicity 
for all three systems, as shown in Fig.l. It is clear that for Au and La 72 > 2 at the peak, 
while for Kr 72 < 2 . The position of the maximum 72 value defines the critical point, m c , 
where the fluctuations in the fragment sizes are the largest. To obtain m c accurately for 
each system a high resolution version of Fig.l, with points corresponding to each value of 
m, was fitted with a polynomial of order 3-9 and the fit with the best chi-square per degree 
of freedom was then chosen to locate the multiplicity at which 72 reaches a maximum. The 
decrease in 72 with decrease in system size observed in Fig.l is also seen in 3d percolation 
studies and has been attributed to finite size effects [0 . 
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FIG. 2. a). Fluctuations in the size of the largest fragment as a function of multiplicity. 

b). Excitation energy as a function of multiplicty. 

Another way of identifying the critical point is from the fluctuations in the size of the 
largest fragment. The fluctuations in this quantity, A^max, peak at the critical point as 
shown in Fig.2(a). For Kr the peak in AAmax is not as well defined as for Au and La. One 
sees a peak in A^ max for Kr at m ~10. This corresponds to ~ 3 MeV/nucleon excitation 
energy, which is too low for MF to occur. Thus in case of Kr the m c value was obtained from 
Fig.l only. For Au and La m c was obtained as the average of the two peak values in Figs.l 
and 2. The m c values for Au, La, and Kr are 28±3, 24±3, and 18±2, respectively. The m c 
value for Au is in agreement with our earlier reported values for Au within the respective 
uncertainities |]J 

The thermal excitation energy, E£ h , i.e. the energy available for particle and fragment 
emission, is a more fundamental quantity than the multiplicity. The experimental relation 
between these two quantities [f| is shown in Fig.2(b). 
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FIG. 3. &).ln{Mz)vsln(M2) for Au above the critical multiplicity, b). The average value of 
ln(Ms) at a given ln{M2). 

To extract the power law exponent r, we examine the MF region above m c , i.e. the 



region past the peak in Fig.l [26]. Fig.3(a) shows a scatter plot of ln(Ms) vs ln(M2) for 
Au [pl,|16|)p2| . The slope, S, of the line through the points is related to the exponent r, 
S—(r — 4)/(r — 3). To obtain r we plot the average value of /n(M 3 ) vs /n(M 2 ), Fig. 3(b). We 
fit the region between E^ h — 5.5 -7.5 MeV/nucleon to obtain the r value. The lower energy 
is ~ 1 MeV/nucleon higher than the energy corresponding to the peak in 72 and the higher 
value is close to the end of 72 branch above m c in Fig.l. We obtain r = 2.16 ± 0.08 for Au 
with x 2 /dof of 1. This value is in agreement with the r value from the single parameter 
fit, ua = %A~ T , at m c [§|,^TJ. The same procedure was followed to fit ln(M 3 ) vs ln(M2) for 
La as shown in Fig.4(b), derived from the scatter plot for La in Fig.4(a). For La we obtain 
t = 2.10 ± 0.06 , with x 2 /dof ~ 6, again in agreement with the value obtained from the 
one-parameter fit at m c . 

The data for Kr are shown in Fig. 5 (a). There is a distinct difference between Fig. 3 (a), 



Fig. 4(a) and Fig.5(a). In the plot for Au, ln{M$) and ln(M2) lie on a very narrow band while 
for Kr there is a large variation. This difference reflects a wider fragment mass distribution 
for different events with the same multiplicity for Kr. A similar trend with system size is seen 



in percolation indicating that this is a finite size effect [ 26| . Fig. 5(b) shows the fit for Kr. A 
nonlinear behavior is clearly observed. This contrasts with the linear behavior seen in the 
corresponding plots for Au and La in Figs. 3(b) and 4(b) respectively. An exponential fit to 
the data is shown in Fig. 5(b) to guide the eye. The fitting region was chosen by the criteria 
laid down in case of Au and La. A linear fit to the Kr data gives a value of r = 1.88 ± 0.08 
with an unacceptably large x 2 /dof of 20. This result is consistent with Fig.l, in which the 
peak 7 2 value is < 2 for Kr. 
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FIG. 4. a).ln{M^)vsln(M2) for La above the critical multiplicity, b). The average value of 
ln{M%) at a given ln{M2). 

The thermal excitation energy per nucleon at m c , E*, was obtained for each system from 
the variation of E^ h with m as shown in Fig.2(b). The dependence of E* on system size is 
shown in Fig. 6 (a), where the size of fragmenting system is the average remnant mass at m c 



|| . The width of the remnant mass distribution at m c is indicated by the horizontal error bars 
and is ~ 6-8% [|,§. Fig. 6(b) shows the isotope freeze-out temperature, T He _ DT , obtained 
from the H 2 /H 3 and He 3 /He A double isotope ratios at m c p,27|. Both E* and Tue-DT 
decrease with increasing system size. We can compare these results with calculations which 
have studied highly excited nuclear matter. The temperature-dependent Hartree-Fock(HF) 
calculations for equilibrated hot nuclei show that Coulomb repulsion causes the compound 
nucleus to become unstable at a lower temperature than the uncharged system p8| . The 
trend seen in the present work is also seen in a HF calculation using a Skyrme interaction with 
a soft equation of state |2£j. This temperature is shown in Fig.6(b) as Tu m i t . In another 



study J3(J it was found that finite size effects and Coulomb force lead to a considerable 
reduction in the "critical" temperature. 
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FIG. 5. (a). ln(Mz)v sln(M2) for Kr above the critical multiplicity, (b). Average ln{M^) as a 

function of ln{M2) ■ 



The Au, La and Kr results can also be compared with the statistical multifragmentation 
model (SMM) [|TQ| , p~T|] . The SMM E* values are shown in Fig. 6 (a). The agreement between 



data and SMM is good although a ~ 1 MeV/A discrepancy is observed for Kr. The SMM 
breakup temperature Tsmm [TT[ is shown in Fig.6(b). There is a decrease in both temper- 
atures with increasing system size. It is apparent that the Tne-DT temperature is about 1 
MeV lower than the SMM temperature. This difference is due to the fact that Tue-DT is 
measured after secondary decay has taken place, while Tsmm corresponds to the breakup 
configuration. Particularly interesting is the fact that the experimental Tu_e-DT tracks Tsmm 
in its dependence on system size at m c . SMM indicates that the decrease in both Tsmm and 
in E* with increasing system size is due to the increase of the Coulomb energy. This result 
suggests that the Coulomb energy plays a central role in the MF of nuclei. 
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FIG. 6. a). Energy (MeV/nucleon) at m = m c . b). Tn e ~DT, Tsmm and Tu m n as a function 

of the system size. The lines through the points are linear fits to the data. 



The micro canonical Metropolis Monte Carlo (MMMC) 0,0] calculations have empha- 
sized that MF is controlled by the competition between long range Coulomb forces and finite 
size effects. Finite size effects in models with only short range forces predict an increase in 
the critical temperature as the system size increases, as is evident from percolation |31] and 



Ising model studies f32| . Since the experimental temperature exhibits the opposite depen- 
dence on system size, it is apparent that Coulomb effects are more important than finite size 
effects. For finite neutral matter the critical temperature (T c ) is expected to be ~ 15-20 MeV 



30| , |33|| . The observed T c for A=160 is ~ 6 MeV. Compared to finite uncharged nuclei, the 



presence of Coulomb energy plays a role in lowering the excitation energy needed to reach 
the regime where critical signatures are observed. In the smaller Kr system there is less 
Coulomb energy in the initial remnant state. Achieving multifragmentation in this system 
requires greater excitation energy /nucleon compared to Au and La (as shown in Fig. 6) and 
as a result, the dynamics of the ensuing disassembly may not take this system near its critical 
regime. 

In conclusion, we have analyzed the fragment distributions resulting from 1A GeV Au, 
La, and Kr on carbon. The reduced variance 72 has a peak at the multiplicity where the 
fluctuations in A max are largest. The peak value of 72 is > 2 for Au and La and they exhibit 
a power law fragment yield distribution at m c . The peak value for Kr is < 2 and this system 
does not exhibit a power law with r > 2. The decrease in 72 with decreasing system size 
can be attributed to finite size effects. These observations argue against a continuous phase 
transition in the MF of Kr but are consistent with such a transition in the MF of La and 
Au. Recent analysis based on the SMM microcanonical caloric curve [|IJ, which indicated 
a first order phase transition for the MF of Kr and a continuous phase transition for the 
MF of Au is consistent with experimental observations. The observed decrease in excitation 
energy and temperature with an increase in system size for MF at the critical point shows 
the importance of the Coulomb energy in MF. 
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